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Abstract 



We show that a low energy beta-beam facility can be used to search for sterile neutrinos by 
measuring the disappearance of electron anti-neutrinos. This channel is particularly sensitive since 
it allows to use inverse beta decay as detection reaction; thus it is free from hadronic uncertainties, 
provided the neutrino energy is below the pion production threshold. This corresponds to a choice 
of the Lorentz 7 ~ 30 for the 6 He parent ion. Moreover, a disappearance measurement allows the 
constraint of sterile neutrino properties independently of any CP violating effects. A moderate 
detector size of a few 100 tons and ion production rates of ~ 2 • 10 13 s _1 are sufficient to constrain 
mixing angles as small as sin 29 = 10" 2 at 99% confidence level. 
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I. INTRODUCTION 



A large number of experiments have now convincingly demonstrated that active neutrinos, 
i.e. left handed neutrinos which interact via W and Z exchange, can change their flavor. 
More recently, KamLAND [[I] is providing direct evidence for neutrino oscillations. The 
oscillation of three 1 active neutrinos describes the global neutrino data very well, see e.g. [[|. 
The fact that neutrinos oscillate implies that at least two of the mass eigenstates have a 
non-zero mass. Most models to accommodate massive neutrinos introduce right handed 
neutrinos, i.e. states which do not couple to W and Z bosons. These right handed neutrinos 
can either directly provide a Dirac mass term or they mediate the seesaw mechanism In 
the latter case, the right handed neutrino tends to be very heavy with ~ 10 12 — 10 15 GeV. 
However, in the most general scenario there is a 6x6 mass matrix whose entries are essentially 
unknown. To obtain the physical neutrino masses at scales below the electroweak phase 
transition this mass matrix has to be diagonalized and its eigenvalues are the neutrino 
masses relevant for low energy observations. Since the entries of the mass matrix are not 
known it can gives rise to any spectrum of eigenvalues and thus neutrino masses. The need 
to describe the oscillation of three active neutrinos implies that at least three eigenvalues 
have to be of (9(1 eV) or less. However, there can be — 3 additional small eigenvalues 
corresponding to light neutrino states, which due to the Z decay width bound would have 
to be sterile. Thus, we conclude that sterile neutrinos are theoretically well motivated by 
the observation of neutrino mass. Furthermore, we see that there can be one or more sterile 
neutrinos which are light enough to play a role in neutrino oscillations. Note, that all 
these considerations are entirely independent of any experimental claims to have seen sterile 
neutrinos, like, for example, the one by LSND ||. 

Interpretations of the LSND claim in terms of only sterile neutrino oscillation have 
failed 0. The remaining models for explaining LSND, see e.g. J/], are viable explana- 
tions of the data only because of the fact that direct tests of LSND, like MiniBooNE ||, 
employ a different baseline or energy, while the ratio L/E is close to the original experiment. 
Therefore, the status of LSND can only be settled by experiments using the same L and the 
same E. 

Light sterile neutrinos have a large number of phenomenological consequences. Most 
notably, they would contribute to the energy density of the Universe and being highly 
relativistic, they would act as hot dark matter. This allows us to use various cosmological 
data sets to put severe bounds on the mass of the light sterile neutrinos of order 1 eV M . 
Among other things, these bounds rely on the sterile neutrinos being in equilibrium with 
the surrounding thermal bath. Since they are sterile, the only interaction 2 they have is via 
mixing with active neutrinos. Thus the cosmological mass bounds only apply if the mixing 
with active neutrinos is sufficiently large. Even very weakly mixing sterile neutrinos can 
have visible effects in astrophysics fl0|-|T3f or Big Bang nucleosynthesis [14, [T5]. However, 



these small mixings can be beyond the reach of terrestrial experiments. 

In this work we explore how well sterile neutrinos can be constrained by a dedicated 



1 We know from the invisible decay width of the Z boson 0, that there are 3 active neutrinos with m < 
m z /2. 

2 besides gravity, which however is far too small for this purpose. 
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oscillation experiment based on a low energy beta-beam facility. For other physics which 
can be studied using a low energy beta-beam, e.g. see |[5J. The proposed experiment is a 



disappearance experiment and will be sensitive to Am 2 = 0.5 — 50 eV 2 and can probe mixing 
angles as small as sin 2 29 = 10~ 3 — 10~ 2 . The oscillation probability for our purposes is a 
two flavor v e — > v s oscillation and the survival probability of v e is given by 

„ 9 Am 2 L , , 

P- e - e = 1- sin 2 2^ sin 2 — — . (1) 

4.C/ 

Note, that this parametrization is phenomenologically complete, in the sense that 
the electon neutrino disappearance probability in any 3 + iV neutrino scenario can be 
parametrized like given in equation [1]. For instance, in the common 3 + 1 scheme, sin 2 29 = 
4|[/ e4 | 2 (l — | t/ e4 |) 2 , for small sin 2 29 this can inverted and we find |?7 e4 | 2 ~ 1/4 sin 2 29. On the 
other hand, sin 2 2# SB L = 4 1 t/ e 4 1 2 1 1 2 , where 9$bl is the angle constrained by short baseline 
appearance experiments like LSND or MiniBooNE. |f^4| 2 is constrained by disappear- 
ance experiments like CDHS and found to be smaller than 0.1 and thus, a bound on sin 2 29 



becomes bound on sin 2 2#sbl < 0.1 sin 2 29, see e.g. reference [17|. 

The most sensitive experiments in this mass range looking for the disappearance of P e 
have been the reactor neutrino experiments: Bugey [IS| and Chooz ]T9| . The Bugey bound 



is valid in the range Am 2 = 0.1 - leV 2 and it can constrain sin 2 29 below 5-10 2 . The 
bound from Chooz is sin 2 29 ~ 10 _1 for Am 2 > 0.01 eV 2 . 



II. EXPERIMENTAL SETUP 



The concept we propose here is based on using a pure i> e beam from the beta decay 
of completely ionized radioactive ions circulating inside a storage ring. The Lorentz boost 
7 of these ions will be chosen such that the resulting v e have an energy below the pion 
production threshold. In this case, the by far most likely reaction is inverse beta decay 
on free protons. In this experiment we have an accurate theoretical understanding of the 
neutrino flux, spectrum and cross section. The detector will be placed so close to the 
neutrino source that oscillation will happen within the detector itself. Thus the different 
parts of the detector will effectively act as near and far detector. This allows the cancellation 
of most systematical errors in a similar fashion as in modern reactor neutrino experiments. 
A concept exploiting several oscillation maxima inside one detector has been proposed using 
very low energy neutrinos from a stationary (7 = 1) radioactive source inside the LENS 
detector. This concept is called LENS-Sterile and more details and sensitivity estimates can 
be found in [p0| . 



Beta-beam 



In [21], the idea of using beta decay of exotic nuclei to produce well defined neutrino beams 



was introduced. The basic observation is, that if a beta decaying nucleus is moving with a 
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Lorentz 7>1, the isotropic 3 neutrino emission will be collimated into a beam. A beta-beam 
facility consists of four parts: isotope production, ionization, acceleration and storage ring. 
Obviously, only isotopes which have the right lifetime are suitable for a beta-beam; they 
need to be long lived enough to allow sufficient time for beam formation and acceleration. 
On the other hand they need to be short lived enough to produce a reasonable neutrino 



flux. Isotopes with lifetimes around Is turn out to be the most suitable |^T[. The isotope 
we consider here is 6 He, which beta decays with a half life of 0.81 s and has an end-point 
energy Eq = 4.02 MeV. The production of exotic isotopes is a well understood technology 
which has been developed for studies of nuclear physics; for a review, see e.g. |22]. 6 He is 



most efficiently produced by the ISOL method in conjunction with a so called converter: a 
proton beam impinges on a primary spallation target, the converter, which in turn produces 
neutrons. These neutrons then hit a Beryllium target where the 6 He is produced. Helium, 
being a noble gas, easily escapes the Beryllium target. With this approach, production rates, 
p, of 2 • 10 13 ionss -1 are achievable |23[ . 



However, the number of ions which can actually be injected into the storage ring is 
significantly smaller and we will present a simplified estimate of this number. First of all, 
there will be decay losses, which are incurred mostly while the ions are still at very low 7. 
The fraction of ions lost to decays depends on the ratio of the beta lifetime and the time 
required for ionization, beam formation and initial acceleration. We assume that this time 
is dominated by the cycle time, t c , of the accelerator and that the decay losses happen while 
waiting for the accelerator. Secondly, there will be the usual losses from beam handling and 
acceleration, e a . We assume e a = 0.5. The number of ions which can can be accelerated in 
each cycle, n c , is given by 

n c = £ a pT (l - e - ^) . (2) 



Also, during acceleration a number of ions will decay and we assume that the 7 rises linearly 
from injection into the accelerator till to the point where the ions leave the accelerator. The 
number of ions after acceleration, n a , is 

_ t c 

n a = n c 7 t-*- 1 ^ . (3) 

The rate of ions, n.j, which can be injected into the storage ring is rii = n a /t c . Note, 
that this expression has the correct asymptotic behavior: lim^-^Q rij = p. Typical cycle 
times range from around 8 s (CERN PS) down to about 1.8 s (FNAL main injector). Taking 
7 = 30, this yields a range for rij = (0.3 — 3) • 10 12 s _1 . Note, that these numbers are based 
on current technology and no special effort for optimization has been made. Therefore a 
rate of 3 • 10 12 s -1 ions injected into the storage is conservative and will be our default setup 
(corresponding to an accelerator like the FNAL main injector). The issue of luminosity 
and useful ions for beta beams has been studied in much more detail and the resulting 
parameters can be found at f24jl . In this considerably more realistic study an injection rate 



of 1.5 • 10 12 s 1 is found, albeit at 7 = 100. Given, that at lower 7 the losses in acceleration 
are less, our simplistic estimate agrees very well with this number. Furthermore, we assume 



3 Beta decay is naturally isotropic for nuclei with initial spin of and for the other nuclei the decay is 
isotropic for an ensemble of ions with unpolarized nuclear spins. 
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FIG. 1: Schematic of the detector and storage ring. 



that the beta-beam complex operates for 1.6 ■ 10 7 s per calendar year and our experiment 
runs for a total of five years. 

In order to compute the fraction of useful decays, i. e. those decays which happen in the 
straight section of the storage ring, we need to specify the geometry of the storage ring. 
We consider a magnetic field B = 5 T inside the storage ring. Taking 7 = 30 this yields 
magnetic radius p = 56 m for the ^He -1-4 " ion. The useful fraction, /, of decays then is given 
by the ratio of the length of the straight section, S, to the overall circumference of the ring. 

/ = • (4) 

J 2S + 2np K J 

In an ordinary beta-beam experiment the goal is to maximize / by choosing S to be very 
large. This is possible because the baseline, L, is much larger than S, L ^> S. Thus the 
storage ring, despite its large size, can be approximated as a point source. Here, however, 
we will consider the case where S and L are of the same order, S ~ L. If the goal is to 
measure oscillation we need to be have an accurate determination of L/E. Obviously, a 
lower bound on the uncertainty on L/E is given by S itself. Therefore, we would like to 
have S small enough as to maintain the oscillatory signature. On the other hand, S should 
be as large as possible to have the largest possible event rate. As we will discuss later in 
detail, a good choice of S — 10 m. Thus, we obtain / = 2.7%. Note, that due to very short 
straight section, the majority of decays will happen in the arcs; this will lead to considerable 
energy deposition in the magnets, which may exceed the heat load or radiation tolerance of 
the magnets. This problem needs further study, but is outside of the scope of this work. 



B. Event rate calculation 

Figure [I] shows a schematic of the setup we consider. In this scheme we will have a 
cylindrical detector whose symmetry axis is aligned with the straight section of the storage 
ring. The free parameters are (see figure [l]): the length of the straight section (S), the 
distance between the front end of the storage ring and the front end of the cylindrical 
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detector (Y), the radius (R) and the length (D) of the detector. In the following, we will 
derive a general expression for the neutrino event rate as a a function of these free parameters. 

Neglecting the small Coulomb corrections to the beta-spectrum 4 , the lab frame neutrino 
beta-beam flux per unit length of the straight section in units of sr _1 MeV -1 s _1 m -1 emitted 
at an angle 9 with the beam axis is described by 

1 



^ar {E Q) = * (E - E*)E* 2 ^(E - E*)* - ml (5) 

47r mP e j 7(1 — p cos 9) v 

where m e is the electron mass, E is the electron total end-point energy and E* is the rest 
frame energy of the emitted neutrino 5 . / is the phase space factor associated with the beta 
decay of the nucleus. 7 is the Lorentz boost such that E* = jE(l — (3cos9), E being the 
neutrino energy in the lab frame, g = N /S is the number of useful decays per unit time 
per unit length of the straight section. 

To calculate the resulting number of events in a cylindrical detector of radius R and 
length D aligned with the beam axis it is necessary to integrate over the length S of the 
straight section of the storage ring and the volume of the detector. Here we assume that the 
beam is perfectly collinear and has no transverse extension 6 . The un-oscillated event rate 
in a detector placed at a distance Y from the storage ring is given by 



dN r s r D r e ' r E ' 

— = ne dx dl\ d92nsm9 dE <f> Near {E,9) a(E), (6) 
at Jo Jo Jo JE min 



where 



tanfl'(x,l)= R f and E> = ^ (7) 
v ' Y + x + £ 7(1 -p cos 9) 

Note, that the baseline, relevant for oscillations, is L = Y + x + I. Here, n represents 
the number of target nucleons per unit detector volume, e is the detector efficiency which 
is taken to be unity in our calculation. E min denotes the energy threshold for our detection 
method. We work with a threshold of 25 MeV which ensures that our events are well above 



the backgrounds. cr(E) stands for the inverse beta decay reaction cross section [25|, which is 
the predominant reaction channel at the considered energies. Figure |^ shows the resulting 
un-oscillated event rates for S = 10 m, Y = 50 m, R = 3.6 m and D = 28.7m. This 
corresponds to a detector mass of 1 kton. 

C. Statistical Analysis 

For our statistical analysis we use the so called pull approach as used in PBL E7I]. We bin 



our data in L/E into bins of equal width. The y 2 function is defined as 

(8) 



i=l 



4 We checked that these corrections are negligible for our purposes. 

5 Quantities without the V refer to the lab frame. 



Note, that the beam size is of order 10 m, whereas all other length scales are of order ~ 10 m. 
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FIG. 2: The un-oscillated event rate as a function of neutrino energy. The event rate has been calculated 
including all geometrical effects and with a luminosty of 3 ■ 10 12 ionss -1 . The different lines show the result 
for different values of the length of the straight section, S, as indicated by the labels next to each line. 

where n is the total number of L/E bins and 



U{ Am 2 , sin 2 26}, {£,}) = N? ({Am 2 , sin 2 26}) [1 + ir%] , (9) 

where Nf h ({Am 2 , sin 2 26}) is the predicted number of events in the z-th L/E bin for a set 
of oscillation parameters Am 2 and sin 2 26. The quantity n s in Eq. (J|) is the systematic error 
on the normalization of our signal. We have taken ir s = 1% to estimate the performance of 
our standard set up. The quantity £ s is the "pull" due to the systematical error on signal. 
In Eq. (U), Nf x corresponds to the data of the experiment. Since there is no data yet, 
this number of observed signal events in the detector has been computed for the case of 
no neutrino oscillation. We assume that our setup is background free, see also section |11E. 



Even if this assumption were violated we do not expect this to affect our sensitivity since 
the background will not be able to mimick the oscillatory behavior of the signal as shown 
in figure |3|. 



D. Optimization of the geometry 

The goal is to obtain an experimental configuration which has optimal sensitivity to the 
disappearance of v e corresponding to a mass squared difference Am 2 = 1 — 10 eV 2 . For this 
optimization we fixed the detector mass at 1 kton, thus D is entirely determined by R or 
vice versa. We tested values of the detector distance Y = {30, 50, 70, 90} m. We studied 
variations in 7 from 20 — 35. This range was chosen to stay below or close to the pion 
production threshold. We also changed the detector radius from 3.5 — 4.5 m. We found that 
within those options the configuration with 

7 = 30 5= 10 m 
L = 50m D = 28.7m (10) 
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FIG. 3: This figure shows the ratio of oscillated to un-oscillated events as a function of the reconstructed 
L/E. In left panel, Am 2 = 10 eV 2 and in right panel Am 2 = 40 eV 2 . The value of mixing term sin 2 28 = 0.1. 
The red (solid) line includes all geometrical effects and the detector resolution, whereas the blue (dashed) 
line assumes a point source of neutrinos. 



is optimal. It is easy to see that for an average baseline (S/2+Y + D/2) ~ 70 m and a 7 = 30 
we have the first oscillation maximum around Am 2 ~ 2eV 2 . Note, that the parameters in 



equation 10 yield both a baseline and energy close to the ones of LSND. 



Using these numbers, figure |3| shows the resulting ratio of oscillated to un-oscillated event 
rates for two different vales of Am 2 as a function of the reconstructed L/E. The blue line 
assumes that the neutrinos are all generated at in one point in the middle of the straight 
section. The red line fully accounts for all the geometry effects. Clearly, for Am 2 = 10 eV 2 
(left hand panel), several oscillation periods can be resolved. A comparison between the 
amplitudes of those periods allows to cancel systematics to a large extent (see also, figure [|). 
At Am 2 = 40 eV 2 (right hand panel) only an average suppression can be observed and the 
sensitivity is entirely determined by the achievable systematic errors. 



E. Detector 



The detector we envisage is essentially similar to the MiniBooNE detector [28], however 



with a cylindrical shape. The important background will all be beam based, since beam-off 
backgrounds will be well measured and cosmic ray events can be tagged with high efficiency 
as was done in similar near surface experiments |]2"§fl . The beam energy has been selected 
to be below threshold for pion production, therefore there are few channels available for 
neutrino interactions. Only charged current quasi-elastic scattering on carbon and electron 
elastic scattering can mimic the signal of the inverse beta decay primary positron. These 
event types will experience the same disappearance rates due to oscillations, but the neutrino 
energy reconstructed under the inverse beta decay hypothesis would be systematically less 
than the true neutrino energy. At these energies the cross sections for these background 
interactions are smaller by at least an order of magnitude, so any effect they might have on 
the measured oscillation parameters, in particular sin 2 26 would be at most 10%, and they 
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can be accounted for and corrected. 

To further reduce the impact of these beam based backgrounds, the detector should be 
optimized for the detection of inverse beta decay. Typically this is done by tagging the 
primary positron with the free neutron capture in delayed coincidence with a mean lifetime 
of ~ 200 /zs in undoped organic scintillator. However, observing the 2.2 MeV gamma ray 
from neutron capture on hydrogen can be a significant challenge in a detector tuned to 
see events in the 50 to 200 MeV range. Additionally, the long delay time will put several 
beam bunches between the primary and secondary events and will therefore increase the 
probability of false tags. Adding gadolinium to the scintillator would help somewhat by 
increasing the tag energy to 8 MeV and reducing the capture time to ~ 30 fis. Nevertheless, 
even if the neutron tag was highly efficient, at least some of the quasi-elastic events on carbon 
will also have a correlated neutron capture tag. Another possibility is to design the detector 
to be sensitive to the positron direction. We expect the elastic scattering events to be peaked 
in the very forward direction, while the quasi-elastic carbon events will have a much broader 
angular distribution. The angular distribution of the hydrogen inverse beta decay events 
will fall somewhere in between. Sensitivity to the angular distribution can be achieved by 
reducing the scintillation light to a point where Cerenkov light can be distinguished. 

To achieve an energy resolution of 10% or better over the 50 to 200 MeV range should 
be possible with a photo-cathode coverage of 10% if we assume approximately equal parts 
Cerenkov and scintillation light |p8 |. Position and timing resolutions, needed to correlate 



events with beam bunches, should be achievable at the half meter and 10 ns level. With a 
bunch spacing of 100 ns or greater, this resolution would provide sufficient space between 
bunches to demonstrate the rejection of non-beam backgrounds. At this level of detector 
resulution the L/E uncertainty is fully dominated by the unkown production point in the 
straight section of the decay ring. 



III. RESULTS 

Figure (| shows the obtainable sensitivity in the sin 2 26 l -Am 2 plane at 99% CL. Our 
default configuration is shown as red solid line in both panels. The setup proposed here 
improves on the existing limit for Am 2 > 0.2 eV 2 . In the range leV 2 < Am 2 < 10 eV 2 
the improvement is one order of magnitude or better. The left hand panel shows how the 
sensitivity changes with varying the length of the straight section. A longer straight section 
(dash-dotted line) implies a large fraction of useful decays and thus better statistics. At 
the same time the L/E resolution is reduced. As a result the sensitivity extends to smaller 
mixing angles (higher statistics) but at smaller Am 2 (resolution). The analogous arguments 
holds also for a shorter straight section (dashed line), which yields smaller statistics but 
better resolution. This improves sensitivity to Am 2 > 10 eV 2 . Thus the length of the 
straight section can effectively be used to tune the experiment to the desired range of Am 2 
and it can be envisaged to have running periods with different straight section lengths within 
the same setup. The right hand panel shows variations of the systematic error tt s as defined 
in equation || Again our default setup with tc s = 0.01 is shown as red, solid line. At very 
small it s = 0.001 the sensitivity (dash-dotted blue line) becomes essentially independent of 
Am 2 , once the first oscillation maximum can be observed. Note, that this very small value 
of tt s = 0.001 probably is not attainable in a real experiment. For more realistic values of 
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FIG. 4: Exclusion plots of sensitivity to active-sterile oscillations in a near detector low energy beta-beam 
set-up. These are 99% CL limits (1 dof). In left panel, we vary the straight section of the storage ring 
while in the right panel, we show the performance taking various values of expected systematic error of the 
considered setup. The pink (dash-dot-dotted line) line shows the current, combined limit on 1 — P gg from 
Bugey pi and Chooz fT§. 
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FIG. 5: Exclusion plots of sensitivity to active-sterile oscillations in near detector low energy beta-beam 
set-up. These are 99% CL limits (1 dof). In left panel, we choose different values of detector mass in the 
range 0.25 to 2 kton. Right panel depicts the performance of the set-up for different values of total injected 
ions. The pink (dash-dot-dotted) line shows the current, combined limit on 1 — P ss from Bugey Q and 
Chooz 



7i s in the range 0.005 — 0.05 the sensitivity limit does not change for Am 2 < 10 eV 2 , which 
is due to the cancellation of the normalization error between different oscillation maxima. 
This is also illustrated in figure |3|. 

Figure [5] shows the obtainable sensitivity in the sin 2 29- Am 2 plane at 99% CL. Our default 
configuration is shown as red solid line in both panels. The left hand panel shows how a 
variation of detector mass by a factor of 10 changes the sensitivity. Remarkably, this change 
is quite small: the detector density and aspect ratio are fixed. Thus with increasing detector 
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mass, the additional detector mass will be exposed to a weaker neutrino flux for purely 
geometrical reasons. The right hand panel shows the change of sensitivity for a changing 
beam luminosity and here the effect is quite pronounced as every additional ion contributes 
equally to improve the sensitivity. Note, the our default beam luminosity of 3 • 10 12 ions s -1 
is based on existing technology and thus may be somewhat conservative. 

IV. CONCLUSIONS 

We have studied a near detector setup at a low-7 beta-beam facility for its ability to 
constrain the disappearance of electron anti-neutrinos for mass squared differences Am 2 = 
1 — 10 eV 2 . The key point is, that for a suitably chosen geometry several oscillation maxima 
occur within the same detector and thus a disappearance measurement at the sub-percent 
level becomes possible without requiring a stringent control on systematic errors. We focused 
on using a beam from the decay of 6 He, which produces electron anti-neutrinos. This 
allows to use inverse beta decay as detection reaction and we can exploit the well defined 
relationship between the positron and neutrino energy. Thus, we have a very clean sample 
of electron anti-neutrino events. We carefully optimized the geometry and beam energy and 
found that 7 = 30 yields the best sensitivity while still having the bulk of neutrinos below 
the pion production threshold. In order to have sufficient resolution in L/E we had to reduce 
the length of the straight section down to 10 m, which makes this setup unique. Note, that 
this allows our experiment to run parasitically in a low energy beta-beam facility since we 
use only around 3% of all ions. Such a low energy beta-beam facility has been discussed 
extensively in the context of using neutrino nucleon scatterng for studies of nuclear structure, 
see e.g. jlSf . For a conservative beam luminosity of 3 • 10 12 ionss _1 and detector mass of 
1 kton we obtain a sensitivity to sin 2 29 ~ 10~ 2 (99% CL) for Am 2 = 1 - 10 eV 2 . 
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